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We study s=| Heisenberg spin ladder with the four spin exchange. Combining numerical results 
with the conformal field theory(CFT) [jjjl, we find a phase transition with central charge c=|. Since 
CN| , this system has an SU(2) symmetry, we can conclude that this critical theory is described by k=2 
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SU(2) Wess-Zumino-Witten model with Z2 symmetry breaking. 
PACS numbers: 75.10.jm 



(N ! I- INTRODUCTION 

Quantum spin ladder systems j|] have been studied in relation with high-T c superconductivity and Haldane's 
O ' conjecture ||. These studies have been done mainly with 2-body interaction terms. Recently from experiments on 
dispersion curve Q|, ||, it is suggested that 4-spin exchange interactions play an important role. 

In general, spin interactions originate from electron exchange interactions. It is well known that Heisenberg model 
, 1 * ' is derived from the second order perturbation expansion in strong coupling limit of Hubbard model. Similarly higher 
order terms give many body spin interactions. For example, fourth order terms give the four body spin interactions, 

e.g. (^Si ■ Sj^j (^Sk ■ Si). In fact, in the one dimensional (ID) two leg ladder, many body spin interactions appear 

from electron exchanges || , jjj . Note that ID nearest neighbor Hubbard model does not give many body interactions, 
because the cyclic path is impossible on this lattice. Therefore, it is important to understand the spin ladder systems 
1 ' including 4-spin exchange term. 

Apart from the spin ladder systems, many-body interactions appear in 2D, 3D systems. It is experimentally known 
that many-body interaction can not be neglected especially in 3 He on graphite ||, and Wigner crystal Q. As another 
example, in high T c superconductor described by 2D Hubbard model, many-body spin interactions appear and they 
may be important physically ]1C[ |. 

This paper is organized as follows. In the next section, we shortly review for s = 1/2 two-leg spin ladder systems. 
In Sec. 3, we discuss the phase transition and the universality class from numerical results and the conformal field 
theory. In Sec. 4, we summarize results. In appendix, we briefly explain theories used in this study, that is, s = 1 
solvable models, WZW model, logarithmic corrections. 



II. OVERVIEW OF SPIN LADDER SYSTEMS 



We treat the following Hamiltonian 

L 



H — Jleg ^ ^ S a ,i ' S at i+i + Jrung ^ <Sl,j ' l?2,j + Jring ^ (Pi,i+1 + -P^i+l) ? (1) 
i—l a =l,2 i=l i=l 

where Ji eg is the intrachain coupling, J rU ng is the interchain coupling, and Pij,k,i is the four spin exchange term which 
is written with spin operators 

P(i, i + l) + p- 1 (i,t + l)= — + Si, v S hi+1 + + S lt i-S 2 , i+ i 
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X. 

H +S2.i'Sl.i+l + Si i-S2A + Si ,i+l-52,i+l 



+4{(S M -S 2 , i )(S M + r S 2 . i +i) 
+(£i,r<S f M+i)(S f 2,rS'2,i+i) 

— (Sl,i-S2A+l){S2,i'SlA+l)}, (2) 

and indexes i are for longitudinal direction, ct=l,2 are indexes of spin chains, (see Fig.0). In this section, we shortly 
review the model and related studies. 
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Quantum spin ladder systems have been well studied. Recently several quantum spin ladder materials are found 
out experimentally, and they attract much attention fs|, [{Ll| . Theoretically, they have been studied with bosoniza- 
tion method and numerical calculations. For only 2-body interaction, there are many studies. See for example, 
§,@, H,[0,i3, ^ references therein. 

When there is the four spin exchange, it is expected that strange phase transitions may appear. In this region, 
relatively there are few studies. Brehmer et al. numerically calculated the dispersion curve in the s=^ spin ladder 
with the four spin exchange ]l6| . They discussed it from the view point of the symmetry breaking. Sakai et al. fTij ] 
studied with numerical calculation for the system with the magnetization plateau. In that paper, they found out that 
this system is Tomonaga-Luttinger (TL) liquid with central charge c = 1, when magnetization m=i, and that BKT 
transition occurs between TL phase and plateau phase region. Honda et al. p7| investigated the energy gap and the 
spin-spin correlation function with the density matrix renormalization group method. Recently experimentists found 
the real material with 4-spin cyclic exchange, i.e., La6CasCu2404i. In that experiment, the dispersion relation was 
observed by neutron scattering ,which suggests four-body terms || , || . 

When J rung < 0, s=tj quantum spin two-leg ladder systems can be mapped to s=l systems. Thus they can be 
considered as a generalization of Haldane's conjecture j| for ladder systems. Haldane phase is characterized by the 
string order parameter, which is related to a hidden Z^ x Z^ symmetry breaking. In s=h ladder systems with four 
spin exchanges, the string order parameter was discussed by Fath et al |18| . In that paper, there are two type gapped 
phases. One is a rung singlet type, another is the AKLT type |l9|. Using density matrix renormalization group, 
Legeza et al. [|15| p ictured phase diagram in these systems. But they have not decided phase boundary properly. 
Nersesyan et al. |2(| referred to the relation between s=| ladder system with four spin exchange and s=l quantum 
spin chain with bilinear- biquadratic term (see Appendix A). In that paper, they examined behaviors of a correlation 
function of relative staggered magnetization and a correlation function of relative dimerization field, and suggested 
the phase transition with central charge c= |. 

When Jmng < 0, we can relate s=^ quantum spin ladder system with four spin exchange to s=l quantum spin chain 
with bilinear-biquadratic term. When J rung > 0, is it obvious ? The discussion of Legeza et al. UJQ is different from 
that of Nersesyan et al. |2^] . Legeza et al. [jl5| wrote that this phase transition is between a nondegenerate singlet 
state and a fourfold-degenerate ground state. The fourfold-degenerate ground state realizes because the singlet-triplet 
gap disappear. Nersesyan et al. [^C| wrote that, there is the region where the triplet Majorana excitation is lower 
than the singlet one, in which it is enough to consider only the triplet Majorana excitation, thus the phase transition 
belongs to the second order type with c = |. Although both of them predicted Takhtajan-Babujian type transition 
point (c = 3/2) in the s = 1/2 spin ladder (with four-body terms), there is a difference on other critical phenomena 
between them. Nersesyan et al. |2(J referred to c = 1/2 (the singlet Majorana field becomes massless in another 
region), Legeza et al. |fj"5| referred to c = 2 Uimin-Lai-Sutherland model |3^| . There remains controversial in this 
region. 

III. RESULTS AND DISCUSSION 

We study Jieg — Jrung— : 1 case with varying J r i ng coupling in the system (1) with no magnetizations, combining 
numerical calculation and conformal field theory (CFT). 

A. numerical results and universality discussion 

At first, we study the size dependence of the ground state energy, in order to see the type of periodicity {Zi or 
Zz) under periodic boundary condition (PBC). In Fig.||, we can see that the ground state energy per site E g (L)/L 
oscillates according to L even or odd (even cases are lower) . This suggests that translational invariance by one site may 
be spontaneously broken. For another evidence of Zi symmetry breaking, we study the dispersion curve, which shows 
that the lowest excitation has soft modes at q — it (see Fig.|J). Thus we will study the system size L=N/2=6,8,10,12. 

In order to discuss what type of phase transitions may occur, we use the effective central charge c, obtained from 
the finite size scaling of the ground state energy. From the conformal field theory at the critical point, the finite size 
scaling of ground state energy under the periodic boundary condition (PBC) is the following |2l|] , p2[ . 

E g (L) = eL-—, (3) 

where E g (L) is the ground state energy for size L, e is the energy per site in the infinite limit, v is the spin wave 
velocity. 
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In addition to the ground state energy, we should obtain v numerically to decide c in equation (3). The spin wave 
velocity v is obtained from 



where q is a wave number. Then we extrapotale v(L) as 



E(q = j-) -£Uy = l» 



(4) 



u (L) = u + + b ( ) + higher order. (5) 



In Fig.[|, we show the obtained effective central charge. The effective central charge shows a maximum value c=1.49 
at Jring— 0.192. According to Zamolodchikov's c-theorem ^3|, the effective central charge shows maximum at the 
infrared fixed point. Thus we can conclude that there is a phase transition with the central charge c = |. 

Only from the central charge, there remain several possibilities for universality class. But considering the symmetry 
of the system, we can discuss in detail. Combining the overall SU(2) symmetry of Hamiltonian, we can relate central 
charge c = | in conformal field theory to the topological coupling constant (Kac-Moody central charge) k in SU(2) 
Wess-Zumino-Witten model. The relation between k and c is given by 

'"ITS' (6) 

Thus the critical theory c=| with SU(2) symmetry is described as the k — 2 SU(2) Wess-Zumino-Witten model. This 
universality class is characterized with scaling dimensions x = | (parity odd,q = 7r), x = 1 (parity even,q = 0) for the 
primary field (in Kac-Moody algebra), x — 2(parity even, marginal, q = 0) for the (Kac-Moody) descendant field (see 
appendix B). They mean that a second order transition occurs with Z 2 symmetry breaking. And the mass (energy) 
gap is generated near the critical point as AE cx | Jring — Jringl fr° m x = 1 term. Here the mass gap AE is related 
with the correlation length £ as AE cx 

Next we study lowest excitations with wave number q = tt, which consist of s = (singlet) and s = 1 (triplet). 
We can expect that these excitation are related with the scaling dimension (x = |,g = it) in the k — 2 SU(2) WZW 
model. In g eneral, from CFT, scaling dimension can be calculated from the excitation energy for finite size under 



PBC as || 



AEi = Et - E = —n. (7) 



Unfortunately, in the k = 2 SU(2) WZW model, there are marginal corrections, 1/lnL, which are difficult to 
treat. Thus we should remove them with the method in appendix C. After removing logarithmic corrections, scaling 
dimension is almost independent on the system size at critical point Jring = 0.192 (See Fig.[|). Therefore, we can 
conclude that this universality class belongs to the k = 2 SU(2) WZW type. 

Finally, we comment on the effect of the marginal operator to the ground state energy and the velocity. For ground 
state energy, it is proved that, 



7TV 

6i V" ' (InL) 



Eg (L) = eL-— (c +7r — 3 (8) 



where o(^l_), o ( ^j^i ) terms do not exist p5fl . We have neglected the o ( q^j^s J contribution (which is small 

enough) to obtain the effective central charge. For the velocity, it is proved that there is not o (1/lnL) term in 
current-current correlation pq| . Therefore, we have neglected logarithmic correction in equation (5). 

B. symmetry breaking and boundary condition 



After studying the universality class, we proceed to discuss the symmetry breaking and the order parameter. 
Although finite size effect prevents us from drawing a definite conclusion, we can consider several clues relating to the 
ordered state. 



4 



At first, for J r i ng < J^ ng , the triplet excitation is lower than the singlet excitation at q = it, whereas for J r i ng > 
J^ ing , the singlet excitation is inclined to be under the triplet excitation. But this is not clear in terms of the finite 
size effect. Secondary, besides the translational symmetry breaking, it is observed a rung-parity symmetry breaking, 
reversing each spin chain {S\ t i <-> 5*2, i), which Brchmer et al. discussed on the dispersion curve. We confirmed 
their dispersion curves at J r i ng — 0.15 (see Fig.^J) . In summary, in the L — > oo limit under PBC, it is expected 
that, for J r ing > Jring, the ground state (q — 0, S 1 " 1 = 0,P rU ng — even) will be degenerate with the lowest state 
(q = it, Sl ot = 0, P„mg — odd), while for J r i ng < J^ , there is no symmetry breaking. 

Honda et al. studied the same system as ours p7|. Using the density matrix renormalization group method, they 
calculated the spin gap and spin-spin correlation function. They claimed that there occurs a phase transition from a 
massive phase to a massless phase at Jring — 0.3. 

Since they studied the open boundary condition, we think that they saw the edge state. For the s — 1 case, in the 
Haldane phase, the edge state can be degenerate under the open boundary condition [P8| . The four lowest energy 
states are almost degenerate, whose energy gap decays exp (— and there is a Haldane gap between them and the 
continuum energy spectrum. On the finite size effect, the edge state is well defined when it is sufficiently localized, 
that is, the system size is large enough compared with the correlation length. However, since the correlation length 
becomes very long near the critical point, so the finite size effect becomes large. Therefore, considering edge states, it 
is possible to explain why Honda et al. estimated the critical point Jr ing larger than ours. Numerically, edge states can 
be confirmed by investigating the ground state energies under open boundary condition with higher magnetizations. 
Note that edge states can be detected experimentally. 

IV. SUMMARY 

In the present paper, we studied s = 4 quantum spin ladder with four spin exchange, in order to see the phase 
transition type in this system. 

We saw the behavior of the central charge and the scaling dimension, which have logarithmic corrections. At first, 
we calculated the central charge, since its correction is small. The central charge has the maximum value c = 1.49 
at J r i„ 9 =0.192. Removing the logarithmic correction for the scaling dimension, we obtained x w | near J r i„ 9 =0.192. 
In this region, the scaling dimension do not have a significant size dependence. So we can conclude that this phase 
transition belongs to a k — 2 SU(2) WZW type. 
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APPENDIX A: S=l QUANTUM SPIN CHAIN WITH BILINEAR-BIQUADRATIC TERM 

We shortly review the following Hamiltonian. 

H = ( COSd fa ■ S n +l) + Sm9 (§ n ■ S n+ l) j . (Al) 

This model is well researched. The phase diagram of this model is known in numerically |p9|| , pb]] . And this model 
is exactly solved by Bethe Ansatz, at the special point 9=—j,j,0. The critical point 9=—j is Takhtajan-Babujian 
point This is the critical point between dimer phase and Haldane phase. It is known that this critical theory 

is described as c = | CFT or k — 2 SU(2) WZW model |3^ ], which is equivalent to three Majorana fermions. 

For general spin s, Takhtajan-Babujian type solvable models are described as Wess-Zumino-Wittten non-linear a 
model(WZW model) with topological coupling constant k = 2s p3], pi], in terms of the analysis for quantum spin 
systems using non-Abelian bosonization f35fl . 

Another Bethe Ansatz solvable point at 9 = j is equivalent to k = 1 SU(3) WZW model with c = 2. This is 
known as Uimin-Lai-Sutherland model |3(|. Above 9 > j , although there is not SU(3) symmetry, the massless region 
extends with the Z 3 symmetry quasi breaking, and soft modes appear at q — ±^-. Recently from non-Abelian 
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bosonization and renormalization group analysis, it is shown that this is a variant of the BKT transition (but the 
different universality class f37|). 



APPENDIX B: WESS-ZUMINO-WITTEN MODEL (KAC-MOODY ALGEBRA) M,\S4\ 



The minimal conformal field theory with the smallest spectrum containing currents obeying the Kac-Moody algebra 
with central charge c is the Wess-Zumino-Witten non-linear a model, with topological coupling constant k. Thus the 
relation between c and k is 

fcdimG 

C = TT^ (B1) 

where dimG is dimension for an arbitrary representation and is called as the dual Coxeter number. In the present 

case, since the system has SU(2) symmetry, we have dimG = 3, and ^ = 2. 

Now we think of the case with the SU(2) symmetry. In this case, primary fields can be classified according to their 
left and right moving spin. There are operators with sl = sr = 0, |, s. Their scaling dimensions are 

The operators with sl half-odd integer (integer) are odd (even) under translation by one site; thus they correspond 
to states with momentum tt (zero). 

Here we think of the case of the k = 2 SU(2) WZW model. In the case of Sl = Sr = |, which forms \ ®\ = 1 © 0, 
the scaling dimension is x = |. In the case of sl = sr = 1 which forms l®l = O0l©2, the scaling dimension is 
x = l. 

Besides Kac-Moody primary operators, there is a marginal operator for all fc, namely Jl • Jr (which is a primary 
field with respect to the Virasoro algebra). Current operator Jl has the conformal weight (h, h) = (1,0) and Jr has 
the conformal weight (h, h) = (0, 1). Thus the Jl ■ Jr operator has the conformal weight (h, h) = (1, 1) which has the 
scaling dimension x = h + h = 2 and the conformal spin 0, i.e. wave number q = 0. 

Current operators Jl, Jr. themselves have conformal spin ±1, corresponding to the wave number q — ±27r/L, thus 
they are related with the spin wave velocity. 



APPENDIX C: LOGARITHMIC CORRECTION TO SCALING DIMENSION (OR ENERGY GAP) 1 34] 



According to the non-Abelian bosonization, in SU(2) symmetric gapless system, the excitation energies (for q = ix 
states), including logarithmic corrections, are 



2nv I \Sl ■ Sr, 
AE = Ei - E w — la* - X ^ / | . (CI) 

where L is a system size, and v is a spin wave velocity. Here S = Sl + Sr is total spin which is a conserved quantity. 
Note that 

1 / - 1 -« 1 - 

Sl - S R = - [S L + S R ) - -S L - -S R (C2) 

= l'(8 + l)-±8 L (a L + l)-±8 R (8 R + l). (C3) 

In the case of sl = sr = |, since si ® sr— h ® h=0 © 1, it forms a triplet s=l and a singlet s=0. The states 
sl = sr = i correspond to a wave vector q = n. For s=l (triplet) case, we obtain 



On the other hand, for s=0(singlet) case, 



(triplet \S L - S R \ triplet) = j. (C4) 



3 

singlet | Sl • Sr | singlet} = — — . (C5) 
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Therefore we obtain 



and 



2irv / 1 1 

AE( S = 1)=E ! -E « — [ Xi --— I, (C(i) 



AE(s = Q)=E i -E ^^(x i + ^^ I -) (C7i 



thus we can remove logarithmic corrections of energy gap. 



Xi = ^- [3AE(s = 1) + AE(s = 0)] . (C8) 
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FIG. 4: Effective central charge 
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FIG. 6: L=12 (N=24) dispersion curve with rung parity at J r i ng = 0.15. 



